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A System for Measuring 
the Thermal Resistance 
of Semiconductor Diodes 

A fast, automatic system for accurately measuring 
junction-to-case thermal resistance of semiconductor diodes. 

By Norman R. Galassi and Bernard S. Siegal 
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Fig. I. This complete thermal resistance measurement 
system speeds evaluation of diodes. In production, it al- 
lows last checking fur packaging faults. The diode under 
test is placed in the heal sink located in the center of the 
rack second from the bottom. 



HliRMAL RP.SISTANCE OF SEMICONDUCTOR DIODES is 

one of the more important parameters specified on 
a diode data sheet. Knowing the guaranteed thermal 
resistance of a particular type of diode, a designer can 
determine if this diode will survive its own internally 
generated heat when used in the equipment he designs. 

The problem in making a determination, however, is 
that many published specifications are not guaranteed 
and often inaccurate, even if unintentionally. Inaccuracy 
arises because of a lack of standardization of measure- 
ment of thermal resistance, and because inaccuracies 
creep in due to rudimentary methods of measurement. 

Attempts have been made to standardize the measure- 
ment of thermal resistance. The Department of Defense 
and the Joint Electron Device Engineering Council 
(JEDEC) of the Electronic Industries Association (EIA) 
have published material toward this end. One such tech- 
nique, called -Method 4081 of MIL-STD-750A: gives 
only an outline of a thermal resistance test. Another 
publication, the Technical Report on the Measurement 
of Thermal Resistance' of J E DEC describes several 
measurement methods and explains the sources of error. 

However, details to implement these methods are not 
fully described. Different manufacturers arc using differ- 
ent techniques which lead to differences in their pub- 
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listed specifications. Therefore, as a safety factor, de- 
signers often reduce the rating from published specifi- 
cations by 50 to 90^ . Diodes with small thermal time- 
constants (e.g., some step-recovery diodes) are the big- 
gest problem since they heat and cool so rapidly that 
accurate measurements are difficult. To resolve these 
problems. Hewlett-Packard Associates has developed a 
system shown in Fig. I . whose accuracy can be demon- 
strated. In addition, the system speeds measurement time 
from the 3 to 10 minutes now required to about 10 to 
20 seconds. 

For most circuit applications, the total thermal resist- 
ance that a device presents to its surroundings consists 
of several separate thermal resistances as expressed in 
Equation 1.' 

0 la - $„ + 0„ + 0,„ ( I ) 

where 

0,„ is the total thermal resistance of the heat path from 
the diode junction to its surroundings. 

0, r is the thermal resistance of the heat path from the 
diode junction to the diode case, 

0 V . is the thermal resistance of the heat path from the 



diode case to the heat sink (mica washer, etc.). and 
0, a is the thermal resistance of the heat path from the 

heat sink to the surroundings (ambient). 
The total thermal resistance (0,.,) is the quantity that 
the circuit designer needs to know. Generally the semi- 
conductor manufacturer has direct control only over the 
junction-to-casc thermal resistance (0/,), and this is the 
parameter usually specified on the data sheet. 

The thermal resistance of any diode or heat sink is 
defined as the ratio of the steady state temperature differ- 
ence between the particular points of interest on the 
device to the power causing the temperature rise.- This 
is expressed by: 




(2) 



where 

0 h . is the thermal resistance of the particular heat path 
under consideration (junction to case). 

T r is the temperature of the cooler end of the path 
(case), 

Tj is the temperature of the hotter end of the path 
(junction). 

P a is the power applied to cause temperature rise. 
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Fig. 3. A typical plot of junction temperature us a func- 
tion of time. Heating power is applied at T. and removed 
at Tt. The slope of several thermal lime constants, t,. t-, 
t,. are shown. These may vary from about 100 ,,s to It) to 
20 seconds. 




Fig. 4. Typical plots of temperature coefficient K (forward 
voltage at I mA) for two types of step recovery diodes. 
Curves I and 2 are for diodes of the same type from dif- 
ferent lots and show thai K is essentially the same, hut 
forward voltage can differ. Curve 3 is for a different type 
of diode, and both K and forward voltage differ from 
Cun'es 1 and 2. 



Measuring Thermal Resistance 

In order to determine 0„- for diodes, it is clear from 
Equation 2 that three quantities must be measured: 
T e , the case temperature, 
Tj, the junction temperature. 

P.„ the power applied to the diode under test in order 
to raise the junction temperature. 

To measure the junction temperature, it is practical 
to monitor a temperature-sensitive parameter. In semi- 
conductor diodes, the forward voltage drop at a small 
fixed forward current (0.01 to 5,0 mA) is an excellent 
parameter to monitor since it is fairly linear with tem- 
perature and is repeatable. 

If the diode under test is mounted in a fairly good 
'infinite' heat sink (0 r , -f- 0, lt small). Fig. 2. the junction 
temperature before the heating power is applied will be 
the same as the case temperature. The case temperature 
changes little during junction heating. 

It is extremely important to remember that the steady 
state temperature difference is used in calculating thermal 
resistance. Most diodes exhibit several separately identi- 
fiable thermal time constants. Fig. 3.' This requires that 
the power be applied for a time much longer than the 
largest time constant of the diode so that the steady state 
condition is reached before the junction temperature is 
measured. In addition, it is necessary to interrupt the 
heating power to measure forward voltage at the forward 
current used to determine junction temperature. This 
must be done quickly after interrupting the heating 
power; the time interval must be small with respect to 
the shortest thermal time constant of the diode under test. 

To measure junction temperature by using the forward 
voltage of the diode under test we must first measure the 
temperature coefficient of forward voltage of the diode 
under test. The temperature coefficient is: 



7". - 7, 



/ml 



I lm = constant 



(3) 



where K is the temperature coefficient (°C/mV) 
7", is the lower test temperature (usually 25°C). 
T : is the higher test temperature (greater than the 
highest that could be reached during the measure- 
ment of $, e ), 
Vfmi is the forward voltage at If m T,. 
Vfma is the forward voltage at /.„, 7\. and 
//„ is the fixed forward measurement current. 
While there is a fairly large (^ 109c ) variation in the 
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(a) (b) 



Fig. 5. Temperature coefficient measurement system (a) is used for characterizing new 
diode types or spot checking during production ij desired. The system (b) (pictured in 
Fig, 1) determines thermal resistance in about 10 to 20 seconds per diode, which enables 
fast measurement o) large quantities of diodes. 



absolute value of the forward voltage at a fixed tempera- 
ture and forward current for diodes from different lots 
(batches), the temperature coefficient is relatively con- 
stant (±3%) from lot to lot for diodes of the same type 
(number). Fig. 4. This allows us to measure the tempera- 
ture coefficient of a sample of diodes of a given type and 
use the average value of the temperature coefficient in 
the thermal resistance calculations for that type diode. 

To heat a semiconductor diode, a known fixed forward 
current is applied to the diode and the resultant forward 
voltage measured. The dc power thus applied is: 

P« = la x Vn (4) 

where 

P„ is the heating power applied to the diode, 
l lh is the heating current applied to the diode. 
Vn, is the forward voltage of the diode at / ,,. 
The measurement of the thermal resistance of diodes 
then requires three steps. 

1. Determine the temperature coefficient of forward 
voltage (K) at the measurement current (/,„,) over 
the specified temperature range (7", and T : speci- 
fied). 

2. Measure the required parameters — 
and Vn at a predetermined / ( ,„ and /.,,. 



where 

V, m i represents the forward voltage corresponding to 
the temperature of the junction heated and 

K,„„ represents the forward voltage corresponding to 
the temperature of the case. 

3. Calculate the thermal resistance by: 



_ (V fmi -V, m .)K 
°" ~ Vn X l n 



(5) 



V. nh V 



Imr 



Thermal Resistance Measurement System 

As shown, the measurement of diode thermal resist- 
ance consists of measuring several different quantities, 
then calculating the thermal resistance of the diode under 
test. This can be a very time consuming task, taking as 
long as three to ten minutes per diode. If an adequate 
characterization of the temperature coefficient is first 
made. Fig. 5(a). the system described in this article can 
calculate the thermal resistance of an average diode in 
between 10 and 20 seconds. This makes the thermal 
resistance calculation system quite useful for evaluating 
packaging techniques, package uniformity, and product 
quality in R & D. production, and Quality Assurance. 
This system has the following pertinent specifications: 
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Fig. 6. The thermal resistance tester portion «/ the measurement system contains the diode 
under test. Programming and timing of the power supply, left, is controlled by the thermal 
resistance calculator and sequencer section. 



■ Heating current (/„,) capability: 0.05 to 5.0 amps 
(Power capability to 40 watts — dependent on the 
diode 

■ Measurement current (/,,„): 1.0 mA 
• Measurement duty cycle: 1.0% 

■ Accuracy: +8. —10% 

The HPA thermal resistance measurement system con- 
sists of three main subsystems: 

■ The Temperature Coefficient Measurement System 

■ The Thermal Resistance Tester (measurement 
circuits) 

" The Thermal Resistance Calculator 
Temperature Coefficient Measurement System 

In practice samples of diodes can be measured on a 
periodic basis to insure that the proper value of K (as 
defined in Eq. 3) is used in measuring the thermal 
resistance. 

The temperature coefficient measurement system used 
for this measurement. Fig. 5(b). consists of an oven, a 
block of large thermal mass (on which the diodes are 
mounted), a constant current source for supplying / lm , 
a temperature measurement system, and a forward volt- 
age readout, such as an HP 3440/3443 DVM. The 
diodes to be measured are mounted in the block and 
placed into the oven. The oven temperature is adjusted 
to the specified 7",. and the forward voltage V, ml is meas- 



ured on all of the diodes at /,„, and 7",. The oven tem- 
perature is then adjusted to the specified ft, and the 
forward voltage V lm2 is measured on all of the diodes at 
/(m and T... The temperature coefficients are then deter- 
mined using Eq. 3. and the average value is calculated 
for a number of sample diodes. 

Since we have found the temperature coefficient K 
to be fairly constant over a wide range of temperatures, 
it might be more practical to use constant temperature 
water baths at 0°C and I00°C to make temperature co- 
efficient measurements on large quantities of diodes. 

Thermal Resistance Tester 

The thermal resistance tester. Fig. 6. performs the 
measurement of V. m .. V, ni . and V,,,. V,„„ and V , m , are 
measured during a 10 /<s interval when the heating cur- 
ent is interrupted and only h m flows through the diode 
under test. V„, is measured during the time the heating 
current is applied to the diode under test. These functions 
arc performed by four main sections: 

1. Shunt switch. The shunt switch. Fig. 7(a), removes 
the heating current /,» from the diode under test 
and the blocking diodes isolate it so only the 1.0 
mA measurement current /,,„ flows. Thus. V and 
V, mi may be measured during the 10 fis interval. 
The shunt switch turns off /,,, for a 10 r s interval 
at a 1,0 kHz rate (1.0% duty cycle). 
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2. Sample gales. The sample gates. Fig. 7(b). are used 
to measure V, mi . V lmt , and V, k . K, m , and V, m , are 
measured by Sampling Gate No. 1. but at different 
points in time (determined by the Thermal Resist- 
ance C alculator). Sampling Gate No. 2 is used to 
measure V t t,. 

3 . Timing circuits. The timing circuits synchronize the 
sampling gates with the shunt switch so they switch 
on and off at the proper times. Fig. 8. 

4. Heating current source. The heating current source 
supplies the heating current la, to the diode under 
test. The heating current source co isists of an HP 
Model 6266 Power Supply programmed to supply 
a predetermined current at a time determined by 
the Thermal Resistance Calculator. 

Thermal Resistance Calculator 

The Thermal Resistance Calculator. Fig. 9. acts as the 
main control unit for the Thermal Resistance Tester, 
current source, and DVM in addition to performing the 
operations required to solve Eq. 5. It consists of seven 
subsections: 



1 . Two sample and hold circuits. One circuit samples 
V .,„. (before the current source is turned on) and 
stores the value in analog form. The second sample 
and hold circuit samples V ',,„, (after variable delay 
interval after the current source is turned on so 
steady state T, is reached) and stores the value in 
analog form. 

2. The subtraction circuit. Performs the operation 
V, mi on the outputs from the sample and 
hold circuits. 

3. The power calculator. Performs the operation V ,„ 
times the programmed /.„, to obtain P„. The value 
of P„ is converted into digital form by the HP 
Model 3440/3443. 

4. The division circuit. Performs the operation: 

0,, K, m , - V lm , 
K ' P„ 

This is accomplished using a multiplying digital- 
to-analog converter and a comparator to divide 
(Yim> — V/mc) in analog form by P„ in digital form. 



From 

HP 6266A 
Power Supr, 



To Timing 
Circuit 




Fig. 7. The shunt switch ( a) of the thermal resistance tester section removes the heating 
current Ir» from the diode under test. The blocking diodes isolate the iliode under test from 
the shunt switch so only the measurement current flows through it. Transistors Q, through 
Q. should be capable of switching 1-2 amperes each in less than 0.5 pS. and the diodes D, 
through D. should have a minority carrier lifetime of less than 0.2 /is at 5 amperes. 

During the measurement interval, the sampling capacitors in the sampling gales ( h) are 
connected to the diode under lest through the FET's. which have a low 'on' resistance. The 
capacitors retain the value oj the voltage sampled since the FETs have a high of)' resist- 
ance at other than the measurement interval, and the high-impedance buffer amplifiers 
isolate the measurement circuitry from the capacitors. 
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Fig. 8. Timing diagram of the 
thermal resistance measurement. 
Sampling Gale #/ samples V (m 
and it can be delayed us long as 
6 /is from the turn-on time of the 
shunt switch. This delay permits 
diodes with a long minority car- 
rier lifetime adequate time to 
reach Vr„, from Vn. Sampling 
Gate #2 samples Vn,. 
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5. The liming circuits. These circuits time all other 
circuits in the Thermal Resistance Calculator as 
well as the power supply and the DVM. Table I. 

6. The programming circuits. These circuits program 
the power calculator as well as the power supply 
and the DVM. 

7. The digital readout. This consists of an HP Model 
3440/3443 DVM. The DVM first reads the value 
of the power input to the diode under test, converts 
it into digital form, and returns this data to the 
Thermal Resistance Calculator. The DVM then 



TABLE I Time Sequence of Operation of 
Thermal Resistance Calculator 



Time 


Action 


To 


Initiate Test 


T, 


Sample and hold V.- c 


T, 


Apply In to the diode under test 


Tto 


Time delay' 


Ti 


Sample and hold V,-,,- calculate 
v\„, - V, mc ; calculate power 


T. 


Read power on DVM 


T, 


Calculate b = V "' ~ V "~ • 
K P. 

Turn ott I*. 


t» 


Read » on DVM 



• The time delay between T, and T, can be set between 1 and 
X seconds to allow the diode under test to reach thermal equi- 
librium before the thermal resistance is calculated. 
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reads the thermal resistance with a scale factor of 
1.0 mV = 1.0°C/watt. 

Errors 

The main sources of measurement error are: 

■ Errors due to variations in the temperature coeffi- 
cient K (estimated to be ±3%). 

■ Errors due to thermal time constants caused by: 

a. Errors in measuring T, before steady state is 
reached. These can be reduced by adjusting the 
variable time delay TD. 

b. Errors due to junction cooling during the meas- 
urement interval. These could be as large as 
— 5% for a device with very short constants 
(100 ,us) and close to zero for a device with long 
time constants. 

■ Errors due to inaccuracies in the Thermal Resist- 
ance Calculator circuits which could be -±2% for 
temperature rises of at least 60°C at the junction 
during testing. 

■ Errors due to T, (case temperature) drift during the 
measurement of T , which could range from 0% 
to +3% for a well designed heat sink. 



Adding these gives a probable overall accuracy of +8, 
— 10%. 
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Fig. L Pushbutton bank On 
tht HP Model 5105 A 
Frequency Synthesizer permits 
selection of high-quality 
signals up to 500 MHz 
in steps as fine as 
0.1 Hz increments. 



Digital Frequency Synthesizer Covering 0.1 

By Alexander Tykulsky 



Frequency synthesizers have become highly versa- 
tile instruments since their introduction as signal 
sources for traditional HF communications. They are 
applied as essential parts of sophisticated systems in such 
areas as spectroscopy, advanced radar, space communi- 
cation and automated testing. 

For some time, there has been a need for high quality 
synthesized microwave signals capable of being switched 
rapidly. Suitable synthesizers, until recently have barely 
reached the VHF range. The use of auxiliary multipliers 
has been the usual expedient for bridging the gap. but 
this is costly and does not lead to optimum performance. 

A new frequency synthesizer, the HP Model 5 1 05 A. 
(Fig. 1) represents a step toward better frequency cover- 
age of the active spectrum. Frequency range is extended 
to 500 MHz, an increase of ten times the range of pre- 
vious HP synthesizers. Design techniques used produce 
these higher output frequencies with less noise than could 
be obtained by multiplying up from lower frequencies. 



The generation of signals in the microwave range is 
also simplified, since by starting at the higher output fre- 
quencies offered by the Model 5 1 05A. less multiplication 
is required and phase noise is appreciably less. 

Synthesis System 

Operation of the synthesis system is accomplished in 
two stages. In the driver unit (HP Model 5 HOB), a 
family of fixed frequencies is derived from a frequency 
standard. In the variable frequency unit (HP Model 
5 105 A), the fixed inputs from the driver are selected as 
required, and further arithmetic operations are per- 
formed to produce the desired frequency. Fixed filters 
are used, and frequency changes are effected only by 
the operation of solid state switches. The method is 
direct! allowing no possibility of an uncontrolled or in- 
correct output. A simplified diagram is shown in Fig. 2. 

The technique of generating a finely-resolved signal 
by means of an iterative sequence of mixer-divider units 
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Ph*u> Modulator 
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Fig. 2. Output from the switched oscillator it processed, carried through two branches and difference 
mixed in the output mixer. Any frequency excursion occurring in the switched oscillator is cancelled out. 



500 MHz in 0.1 Hz Steps 



is borrowed from the HP Model 5100A. 1 The signal so 
generated, covering 30 to 3 1 MHz in decimal steps from 
0.01 Hz to 100 kHz, is added into the N (narrow band) 
branch of the synthesis system at Mixer 1 and is multi- 
plied by 10 in a doubler and quintuples As a result the 
coverage contributed by the mixer-divider units as seen 
at the output is increased to 10 MHz and the steps range 
from 0. 1 Hz lo 1 MHz. A search oscillator is included in 
the mixer-divider sequence to provide continuous tuning 
coverage or FM. 

A signal in the 30 to 39 MHz range, selectable in 
1 MHz steps, is added into the same branch at Mixer 2 
and contributes 10 MHz steps at the output. One strip- 
line oscillator of a set of five (covering 390 to 350 MHz 
in 10 MHz steps) is switched on as required. Simulta- 
neously the spectrum generator is switched to provide 
a complementary signal so that the difference (330 MHz) 
required to activate the N branch is always present. The 

1 Victor E. Van Duzer, '* 0-50 Mc Frequency Synthesizer,' Hewlett- Pitur! Journal, 
Vol. 15. No. 9. May 1964. 



switched oscillator frequency is multiplied by 10 in the 
doubler and quintuplcr of the W (wide band) branch, 
resulting in the 100 MHz step of the synthesizer. 

Finally, a difference is taken in the Final Mixer and 
this difference signal is amplified and leveled. 

Drift Cancelling 

Non-synthesized signals (not quartz crystal derived) 
from the switched oscillators can be used because the 
synthesis system is in a drift cancelling arrangement. 
Cancellation of frequency errors is accomplished by 
channeling the frequencies through two paths as shown 
in Fig. 2. Since any UHF oscillator frequency excursion 
is carried through both paths, it cancels out in the output 
mixer. There is a difference in transmission time through 
the two paths, however, and the correlation is not perfect 
for rapid variations. The potential problem is avoided 
by using oscillators which have no appreciable high fre- 
quency modulation. 
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Fig. 3. Two of the switched oscillators are shown here 
(the two-layer hoard is disassembled). They were designed 
for low noise by placing the active element (a transistor) 
at a relatively low impedance point on the resonator and 
designing the resonators as folded transmission lines in 
a rigid structure to reduce microphonics. In addition, no 
high impedance point is brought out o) the structure and 
oscillator switching is arranged so that resonator Q is not 
degraded. 




Fig. 4. Typical performance in switching from 399 to 
401 MHz. Lower trace is the switching signal. Sweep 
speed is 5 ils per division indicating switching speed is less 
than 20 fa. In photo, the signals are beat with 400 MHz. 
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Switched Oscillators 

Exploitation of the drift cancelling scheme depends on 
the availability of easily switched, high quality oscillators 
(Fig. 3). 

Resonators are tapped transmission lines folded 
within the laminated boards. Tuning is accomplished by 
isolating portions of one ground plane. Output coupling 
results from the proximity of the resonators to the 50 Q 
output strip line. No high impedance point is brought 
outside of the ground planes. The resulting compact 
structure yields the necessary high spectral purity, low 
microphonic susceptibility, and good long-term stability. 

Output Circuit 

Output level is detected at an internal point directly 
on the output transmission line and compared to a refer- 
ence voltage. Control current is developed through a 
temperature compensated dc amplifier and applied to the 
PIN attenuator associated with the final mixer. 

This method of leveling, while permitting efficient use 
of the amplifier output capability, makes the output 
power dependent upon the VSWR introduced by the 
combined effect of the load, transmission line, and con- 
nectors. At the calibrated output level of 0 dBm and with 
an ideal load, the flatness of output over the frequency 
range and for specified environmental conditions is ±: 1 
dBm. Harmonic distortion is typically less than 4%. 

External manual control of the leveling reference pro- 
vides an output capability of approximately —6 to -f 6 
dBm. A metering circuit indicates normal operation of 
the leveling circuit. 

Spurious Outputs 

Non-harmonically related spurious signals are more 
than 70 dB down from the desired signal. Spurious sig- 
nals close in to the desired signal tend to occur in PM 
sideband pairs. Those more than a few MHz away tend 
to be single, additive signals. 

Switching Speed 

Programmed switching may be obtained by use of the 
three remote control connectors in the back panel of 
the variable frequency unit. Switching speed is nominally 
20 /ts. depending upon the frequencies involved. One 
example using the switched oscillators and spectrum 
generator is shown in Fig. 4. 

Phase Modulation 

Phase modulation is accomplished at Mixer 1 by con- 
trolling the tuning of the output circuit. The usable mod- 
ulation rate extends from 0 to more than 1 MHz. with 

'ackard Co. 



the deviation typically down 1 dB at 1 MHz. Variation 
of deviation with output frequency is less than 2 dB. 
Harmonic distortion is a function of modulation depth 
— at 1 radian 5% is typical. Sensitivity is nominally 3 
radians per volt. 

Phase Noise 

Spectral plots of phase noise at 50 and 500 MHz are 
shown in Fig. 5. Also shown are the results to be ob- 
tained by multiplying the HP Model 5100A in noiseless 
multipliers (X 10). It is clear that while the HP Model 
5100A is to be preferred below 50 MHz. the HP Model 
5 105 A is preferable to multiplication (by 10) of the HP 
Model 5100A. The advantage comes from the smaller 
magnification of independent noise originating in the 
filters of the driver unit. This is made possible by the 
use of the drift canceling scheme. 

Solid State Modules 

The design is based upon a large variety of solid state 
circuits, ranging from diode switches to S-band quintu- 
ples. Modular construction is used throughout to achieve 
the required isolation between circuits, to facilitate exten- 
sive production testing, and to simplify servicing. 
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Fig. 5. Comparison of spectral distribution of phase noise 
sidebands of the 50 MHz Model 51Q0A and the 500 Mil; 
Model 5 105 A. The lower frequency synthesizer provides 
improved noise performance at the lower frequencies, 
while at higher frequencies the Model 5105A shows sub- 
stantially belter noise performance than can be obtained 
by multiplying the lower frequency by ten. 



SPECIFICATIONS 

HP Model S10SA Frequency Synthesizer' 
OUTPUT FREQUENCY 0 1 lo BOO MHi 

DIGITAL FREQUENCY SELECTION 0 1 Mi Ihroogf 100 MMf 
p«' Selection trt Tfunl panel pushbutton o' by remote 

a*<un closure An? cnartge in frequency m«r <»• •ccom- 
pliihed in >o i.s typically 

OUTPUT VOLTAGE 

FIXED 0 dBm - l dBm mlo ■ 50-ohm (eaitbve load 
VARIABLE -6 dBm to < 6 dBm into * 10-ohm rwiM 
load typical at 30*C 

SEARCH OflCILLATOR'FM A mrch oscillator provide* con. 
tinuowa variable frequency selection wiln an incremental 
range of 1 0 Hz through 10 Mm/ Manual, or eilernai voll- 
■0* 1*1 *• — M voltal eonlrol arilh linearity of • 5% 
Uanmutn une w«<v rate if ■ (email >* 1000 Hz 

PHASE MOOULATION i50V -oar pane! tOfttfl -3 PJOMI 
meoimum deviation* dc — I MM/ rata Typical aenellivily 
• 9 rad'arr* per volt Linearity depend* woon deviation 
amplitude 

SIGNAL -te-PHASE NOISE RATIO M. a » u .e.1 .n a 30 *Hl band 
'.entered on ma atonal te» eluding a 1 m* band cenfared on 
the tigneli it graaux than 



Output Frequenci — MHz 


1 0D 100 SOD 




41 4* 41 40 



SIGNAL lo-AM NOISE RATIO: (Above 100 iKu GfaaUi man 
74 dB m a 30 »Hz band 



RMS FRACTIONAL FREQUENCY 0EVIATION iwlth a 30 «Hl 
none bandwidth) " 





Output Frequency 


Averaging Tim* 


1 MHi SO MHz 100 MHz 500 MHz 


10 milUaecondt 
t mcond 


1» 10 • ?* 10 ■ t a 10 • 6i 10 
2i 10 • 4 ■ 10 " J • 10" 1 ■ 10 " 



SPURIOUS SIGNALS; Nor. Mrmomcatlv related tignala ara 
at «»aai 70 dB beiow the selected "equencv 

HARMONIC SIGNALS' 25 dB below me aeiected ••eow«nr» 
lapphcabla lo lited output when terminated in 50 ohm*) 

WEIGHT Nel. 75 1|>» |M kgi Snipping. 1U ItH (eO agi 

PRICE I9.1S0.00 

* ttewrn 5111)8 0*i**». 

" tu>i*- HMi tn i!l00 Mm ul.t.jai •■ »"T<ibji fffquarut irAne 
mt, tfua aill iHact IM ilabiUtr tad toKlrai pgnh ol th» 
vvltivl P««l(«m t Q(» Cj< • MM pyj tiM aa ■nta'aal tie 

9M*Cf itaaetrg or tai.cjt« i.-'«fi rr' (oajnpiiliw lo OKf a" 
reHBrtianct a-!» attatail tHM«>« 



HP Model 5110B Synlhesizer Driver 

OUTPUT FREQUENCIES C'i»idti 77 ftftM tieQuancie* lot 
the &ID5A Frequency SyntneaKar 3 0 mrough 3.1 V"r in 
0,1 MHx tlopa 150 mv * t. -> dB) 30 m'oug" 30 MHi tn 
TMHi atecx 24 M»i and 20 MHi H00 mV -IS dB>. 50- 
onm ayitem I MHi buttered output (1 »olt * ' 5 dB Into a 
50-c*i* ie»ia!i»* (oedi erasable ai teat panel conn»ic4o# 

INTERNAL FREOUENCT STANDARD 

TvPf t MHz Quartz OaciHalor 

AO'NO RATE UMM MR - 3 pan* in IC pe» 74 hour* 
STABiliTt . A* a function ot ambient tamperaime 
; 1 1 10 ■ tft from O^C lo 4 »'C 



At a lunction of Ima Miuga • 5 i tO ' tot a • 10% 
change m hne voltage (rated at "5 or 230 vo'ia mi 
line voilagel 

RMS 'PACTIONAL FREQUENCY DEVIATION. 

Average Time 1 MMf Output Frequent,-* 

10 milliaeconda 6 A 10* * 

1 lecond 1 a 10 " 

signal to-PHA5E noise ratio • -as an 

SIONAL-toAU NOISE RATiO ■ >90 08 

HARMONIC SIGNALS ..40 dB bale* the oulpul |*iw 
[■>■ 1 1" Term (net tool 

' IA * SO Mt| t*Ni ( eattWM M tw nnm. •iclat-lf ■ : M; btri 
(tMuM an t*t taiftet 

PHASE LOCKING C*PARIttTY A voltage control feature 
altona 5 oarli m tre treqi^neif confyl lo# -ft lo *S 
■oils applied eitemaii* 

EITERNAL FREQUENC* STANOARD 

INPUT REQUIREMENTS t MHz c 4 MMf 0 7 V rm* Mr* 
5 V maaimum acroat 500 onm« SUbllily and 
apecirai rvuntv o* 5105A Frequency Syntheener will b» 
baniall* determined mr Pv* characteriatica o» Via eolemai 
•tanjufd it uaed 

GENERAL |Af>r>l*« lo SIOSA/SHOBl 

OPERATING TEMP RANGE 0 to * 55 C 

WEiGHT Net 54 lb. |?4 hgi Shipoina. 67 IM (30 Hgi 

POWER "5 o- 230 Will -'0% 50 to 4OT crde* 35 waltf 
eeC* wMtl ht>depe«ideAt *u|K>iie*i 

PRICE M 350 00 

MANUFACTURING DIVISION HP FREQUENCY 4 TIME 
1501 Page Mill Road 
Palo Alio. California 94304 
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Phase Noise in Frequency Synthesizers 



A FREQUENCY SYNTHESIZER OPERATES on the signal 
provided by a frequency standard to produce a de- 
sired frequency which reflects die accuracy and stability 
of the standard. In effect, the frequency of the standard 
is multiplied by a selectable synthesis factor — a rational 
number usually based on the decimal system. The case 
and speed of selecting the synthesis factor, the resolu- 
tion and range of available output frequencies, and the 
spectral purity of the output signal are usually the critical 
characteristics of the instrument. 

Circuits performing the necessary arithmetic opera- 
tions introduce random amplitude and phase variations. 
Amplitude modulation tends to be reduced by subsequent 
non-linear, compressive circuits; phase modulation how- 



ever, remains, and is usually enhanced. Synthesizer noise, 
then, is mainly random phase modulation. 

The total phase noise of a synthesizer includes the 
noise originating in the frequency standard. A useful 
model. Fig. 1(a), shows a simple synthesizer as the com- 
bination of an ideal multiplier and a phase noise modu- 
lator acting on the output of a frequency standard. This 
results in a total modulation made up of two parts w hich 
can be treated separately: The frequency standard noise 
multiplied by the selected synthesis factor, and the sys- 
tem noise <b,(t.F). The system noise results from the 
combined effect of the arithmetic circuits (mixers, mul- 
tipliers, etc.). The importance of a noise source depends 
on its magnitude and the multiplication which follows. 



FREQUENCY 
STANDARD 
Frequency = "<„ 
Phase Noise = putO 




Frequency = Fm» 0 

Phase No.se = F«\,(t> -t-Ml.F) 



(b) 



1 MH.- 

FREQUENCV 

STANDARD 
Frequency = F» 0 
Phase Noise = p„(I) 



390-400 MHl 




y— 0-SOMH7 
Frequency = |F 4 - F B )" t 

Phase Noise = (F. - F,)ffl„(|) +*,(t.F) - « B (!.F) 



390 350 MHz 



FREQUENCY 
STANDARD 
Frequency = ^ 
Phase Noise = <v c (t) 




0-SOMH) 

^ Frequence 

Phase Noise = Fo 0 (t) + 0,(1) 

Nole •>,(!) F) in Fig 1 (a) 
(or the same value of F. 



Fig. I. Synthesizer noise models. Narrow band syslem (a) is a basic form of a typical 
synthesis system. The wideband system (b) is a combination of two narrow-band systems. 
The wideband syslem simplified (e) in which <f>M) = -PxiUF) — <Mt,F) and is virtually 
independent of the synthesis factor F. The numbers given arc examples only. 
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Fig. 2. Distribution of phase noise components. These 
curves also illustrate the sideband distribution if the ordi- 
nate is considered to be two times the single sideband to 
currier ratio, and the abcissa represents offset from carrier. 



Wideband output ranges are most readily obtained 
by down-converting the outputs of two higher frequency 
sub-synthesizers* each of which operates over a prac- 
tical, small percentage bandwidth. The larger synthesis 
factors of each branch lead to greater system noise than 
would be expected from a non-heterodyne technique. 

An example of a wideband synthesizer is shown in 
Fig. Kb). The total phase noise modulation. <i> T (t), to be 
expected is 

(F, - F,,)^(t) + F A <j, ,(t) - F H cb,,(l) 

where the first term is directly related to the frequency 
standard, and the last two terms. constitute systems noise. 
F t and Fa are the selected sub-synthesizer factors, and 
F A — F„ is the overall synthesis factor. For simplicity, 
the model of Fig. Kb) can be represented as shown in 
Fig. Kc) and the system phase noise can be written 

<fr&) = P*tt 4- feft) 

Typical spectral distributions for the model of Fig. 
1(c) are shown in Fig. 2. where <f>.(/). <J>.~(/) and <]>/(/) 
are the Fourier transforms of the corresponding time 
functions. 

Noise derived from the frequency standard becomes 
significant as F is increased, and then predominates at 
very low modulation rates. The best frequency standard 

* SuB-synthesiiets in tnis conteit are the synthesis systems, separate but related, 
container! within a frequency synthesizer 



available, consistent with the cost of the synthesizer, 
clearly should be used. Generally, a crystal filter is em- 
ployed as part of the standard to improve the spectral 
purity. 

System noise is virtually independent of F and is the 
major contributor at rapid modulation rates. Location 
of the corner, ft: (Fig. 2) is dependent upon the narrow- 
est filter bandwidth in the system. In systems where very 
fast switching is not required, the corner fe may be 
moved in by the use of a phase lock loop as a narrow 
band filter, but at the cost of some increase of close-in 
noise. 

The proper design of mixing circuits is a critical 
problem. Low signal levels lead to poor signal-to-noise 
ratios and phase noise. High signal levels may introduce 
spurious modulation at discrete frequencies — an equally 
undesirable result. Mixers applied at the end of a syn- 
thesis sequence introduce less noise since no multipli- 
cation follows. 

Drift cancelling offers another opportunity to reduce 
system noise. The high sub-synthesizer frequencies are 
obtained by adding coherent signals into each branch. 
The non-coherent noise is reduced (since its synthesis 
factor is smaller) while the coherent noise is canceled 
in the output. 

The sub-synthesizer factors for Fig. 3 arc about one- 
sixth of those shown in Fig. 1. The system noise will 
also be about one-sixth — if the noise associated with 
the oscillator cancels. A differential group delay through 
the filters in paths A or B makes the cancellation im- 
perfect. The net effect of the oscillator phase noise can 
be seen at final mixer. Arriving through path A the phase 
is <p„(t + t), where is oscillator noise, / is time and r 
is group delay due to the filters. Through path B the phase 
is </>„(')• The mixer output shows the phase difference: 

<M' + ") - 4>»U) 

This equals -<£;,(/) if r is sufficiently small. Writing dp„(/) 
as the Fourier transform of £„(<). j2-f<f)(f) is the spec- 
trum of <j>:,(t), and j2-xU<$>n(f) is seen at the mixer output. 

It can be seen that if the synthesizer system noise is 
to be minimized $„(/) must be small, particularly for 
large values of /. 

This reasoning also explains the need for a sharp filter 
in the frequency standard output. In Fig. 1(b). assume 
a probable differential delay, r. between sub-synthesizers 
A and B. Then an additional noise modulation tF i,<pM) 
appears in the output and this leads to the additional 
spectral component jTF„2-zf<$)„(f). For sufficiently large 
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390 400 MHe 



FREQUENCY 
STANDARD 
Frequency = F*u 
Phase Noise = « 0 (l) 




0 50MH; 

Frequency = 

(F. - F,)». 

Phase Noise = 
(F« - FeM, (!) + 
*«(t.F) - tttf.n 



Fig. 3. Noise model jor a wideband .synthesizer employing drift cancelling. System noise 
is less than in Fig. 1(b) since F 4 and F« are smaller leading to smaller values of and 
<Pb. The numbers given are examples only. 



values of /, this component becomes overriding and must 
be reduced by filtering. 

Many different synthesis schemes can be devised to 
obtain a given set of numerical results, but the perform- 



ance of each is unique. Choice of a suitable design is a 
compromise involving consideration of discrete spurious 
outputs, switching speed, size, cost and other character- 
istics in addition to phase noise. ■ 
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duplex communications, a pattern generator and a sweep 
oscillator. 

Al designed the UHF section of the Model 5100A Fre- 
quency Synthesizer and was responsible for the design of 
the new Model 5105A. He has developed techniques for 
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Bob Maldewin. right, joined Hewlett-Packard in 1960 and 
participated in the development of the HP Model 5100A 
Frequency Synthesizer. Among other contributions, he is 
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of the Model 5105A stripline oscillators, the stripline out- 
put mixer and the PIN attenuator. 
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Transform Methods for Linear Systems 



By Michael O'Flynn 

Professor of Electrical Engineering, San Jose Stale College 



This is a highly condensed collection of reference material 
on transform methods. Formulas for direct and inverse 
Fourier. Laplace, and Z-transforms are given in compact 
tabular form, with examples of how they are used to find 
system responses for several types of inputs. 

Transforms arc powerful tools for solving linear time- 
invariant system problems, such as 



Given the Input »(f|- 



H(j), the system function 
Of 



lind the output fit) 



the impulse response 
" the system 



Depending upon the system and the excitation, the most 
convenient transforms to use might be Fourier transforms, 
one- or two-sided Laplace transforms, or one- or two-sided 
Z-transforms. These transforms arc treated in many text- 
books (such as references 1 and 3 in Bibliography), and no 
attempt has been made to give an exhaustive treatment here. 
To the author"s knowledge, however, there is nothing in the 
literature equivalent to Table I, which is a comprehensive 
collection of the defining formulas for all of these trans- 
forms and their inverses. The arrangement and notation of 
this table are designed to demonstrate the similarities and 
interrelationships between the different transforms. A fea- 
ture of the treatment of inverses is a logical progression 
from the fundamental inverse formulas to their equivalents 
in terms of complex variables (contour integrals, residues. 
Laurent series). 

Table II applies transform methods to the problem of 
finding the output of a linear system for five types of input 
signals. In each case the transform used is the most appro- 
priate one for the type of input. 

These tables can be used as a reference by anyone work- 
ing in the field of linear systems or — along v. ith a textbook 
— by students of linear system theory. As an aid to the 
student. Table III lists the impedances of the basic network 
elements R, L. and C for the five classes of inputs. 

Classes of Inputs 

In the tables and discussion that follow, five classes of 
inputs are considered. They are: 

i) A damped sinusoidal input x(t) = Ae"J cos (,„„»• + </>„) 

j 

ii) A periodic input x(t) = £ -4„ cos (n,„„t + <•/>„) 

n - 0 

iii) An aperiodic input which has a Fourier or a Laplace 
transform. 

iv) A sampled input .r"(f) known only at specified points. 



v) A random input, where .v(r) is a member of an ergodic 
random process for which the autocorrelation function 
R T h) and its Fourier transform S..(/<„1 are known. 

As an example of the use of the transforms, the aperiodic 
case (iii) will be discussed in detail. The sampled case (iv) 
is similar to the aperiodic case except that Z-transforms are 
used instead of Laplace transforms. The first two cases, i.e., 
periodic and damped sinusoidal inputs, are well known. 

For random inputs (case v) the problem can be restated 
as follows: Given the input autocorrelation function /?,(-) 
and the system function, find the output autocorrelation 
function /?„(-). With this modification, the problem is similar 
to the aperiodic case. 

Inverses by Complex Variable Methods 

The problem of finding the output of a linear system 
whose input and system function are known may be re- 
garded as an exercise in finding inverse transforms. Inverses 
can be found either by looking them up in a table of trans- 
forms, or directly, by computing them using the inverse 
formulas given in Table I. The direct approach requires a 
knowledge of complex variables and contour integration, 
but it gives more insight into the physical meaning of the 
solutions.* For example, the inverse expressions contain 
full information about the natural frequencies of the system. 

Linear System Output for Aperiodic Input 

The output of a linear system for an aperiodic input 
x(t) is 

y(t) — ~ - J ms)X(s)e"ds. 



Here the Laplace transform has been used. /V(.v) is the 
system function and Xts) is the transform of the input x(t). 

If the one-sided Laplace transform is used, and if Jordan"s 
Lemma is satisfied (usually this means simply that 
H(s)X{s)->0 as s-* cc ) then 

y(t) = 0. t < 0 

>'(;) = sum of the residues of the poles of His)X(s)e'', t > 0. 

The complete response y(t) consists of two parts. One is 
the forced response y ,(/). in which only the input frequencies 
appear. [The input x{t) is often called the forcing function.] 
The other part of y(t) is the transient response y,,{t). 

* If s transform has an Infinite number of poles (e.g., a discontinuous function, such 
as _TT_) it may be difficult or impossible to apply the complei-nariable formulas. In 
such cases, time-domain analysis may be simpler than transform methods. 



Reprints of this article are available. Write 
Editor, Hewlett-Packard Journal, 1501 Page Mill 
Road, Palo Alto, California 94304. 
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Table I Definitions of Different 
Transforms with their Inverses 



Fourier Transform 



Laplace Transform (Two-sided) 



Laplace Transform (One-sided) 



Z-Transform (Two-sided) 



/ Hi)e""dl 



Ftim) is defined for ... real and 
exists if 



: J /(/)<- "ill 



/(f) | di exists 



J = <r + j«> 
F(s) exists if 

/I IU)e " ill exists 

This integral exists in some allow- 
able region of convergence 

a, < a <a 3 



S — a -f- /m 



/'(.v) exists if 

/ | /(/)e " | </< exists 



This integral exists in some allow- 
able region of convergence a > a, 



/*(/) = •£l( ll )S(t -nT) 
■i i 

n 

or 

F(z) = F*(s)\_ „ 



This exists in an allowable annulus 
of convergence r, < |z| < r. 



/(f) = 



/(f) 



I /""•••' 



which may be evaluated directly 
or from tables 



Id 



= TT / f-'d^e'-'dd.,,) 



"i < Wo <"-• 

For f < 0 



</5 



Fish" 'tis 



S = /o. 



For f < 0 



/(/) = ^-rif)F(s)e-'fls 
"' C, 

(Sum of residues of 
poles in right half 
plane. 

For I > 0 

l(.l) = -^j(pF(s)e"ds 

= Sum of residues of poles in 
left half plane. 



Sum of residues of\ 
poles with real parts J 
10 right of it,.. / 



For 1 > 0. 

= Sum of residues of poles 
with real parts to left of a„. 

Note: Residues are defined for 
counterclockwise contours. Con- 
tour C, is clockwise. The integral 
around C- is. therefore, minus the 
sum of the residues of the enclosed 
poles. 



/(f) =~ / F(s)e"ds 
For 1 > 0 

/(f) = TZf 'J l'sh-ds 

= Sum of residues of poles of 
F(s).[C encloses all poles of 
F(sl] 

Note: As a consequence of the 
definition of the one-sided trans- 
form an acceptable n., is automat- 
ically to the right of the poles of 

F(s). 



/(«) = -^j'f:Fiz)z"dz 



where C is a circle of radius r„ 
such that r, < r„ < r- : . 

For « ^ 0 

j(n) = sum of residues of F(;);" 
inside C. 

For n < 0 

/(n) = sum of residues of F(:J:." 
outside C. 

Note: the formula for f(n) is the 
formula for the coefficients of the 
Laurent scries for FCt) (see Ap- 
pendix). 





»»f» 0' Conwj»nt€ 







V— 
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Z-Transform (One-sided) 



F(z) = £ /(«); 



F(z) = F*Uh z = e . T 



This exists in an allowable region 
of convergence outside u circle 
centered at the origin. 



10,1 



' /V..-..:- 

2-1 J 



( 



where C encloses all the poles of 

F(z)z" ' 

f(n) = sum of the residues of 
F(z)z n ' 

Note: tile formula for /(/;) is the 
formula for the coefficients Of the 
Laurent series for F{z) (see Ap- 
pendix). 



Im(/I 




Table HI 

Impedance Table 




i ype or inpul /|N 


Transformed Input 


V (51 

Impedance ol Network Element. Z|s) ■ 








R 


L 


C 


i) Damped Sinusoidal 

A t '•' cos(,., f + V I 

Special Cases 

a) Sinusoidal, a, = 0 

b) Exponential. 

„„ = 0. = 0 


At* 

A >■ 
A 


pbasoi 

' notation 


R 
R 
H 


u 

s - — o„ -r /». 

U 
L\ 

'* = — «a 


1 

<« 

l» = — it. + l-« 

| 

1 

0 

i = — a* 


c) dc case. n„ = 0. 

-„ = 0. 4* = 0 


A 




H 


short 
circuit 


1 

°. = o 

open circuit 


ii) Periodic 

V.-f.cos(/w + .(..) 


Same as case i) using superposition 

Each harmonic has an impedance associated with it 


as in (1) (a). 


m) Aperiodic i(f) 


Fourier Transform 

Kji,) = fifi&0 'di 

if integral exists 


R 


!* = /» 


1 




1 .ipl.ice Transform 

/{.r)= J i'(t)e "<li 


R 




1 

Cj 


ivl Sampled f*(«) 


System functions for sampled data systems are usually computed 
first as ratios of Laplace transforms and then converted to 
Z-lransforms. 


v) Random r(f) 

Autocorrelation 
Function — Rjrt 


Power Spectral 
Density 
S,<i.) = 






> i ' 




J l< '->• "di 




N ( «e:S r 0..) = |Z(/-)|'5.(/..) 



Only the poles of .V(.v) contribute to the forced 
response. The poles of His) give the natural fre- 
quencies of the system and contribute to Ihe 
transient response. If we denote the poles of X(s) 

as j,, s.j ,v„ and the poles of 15f(j) as 5,. S., 

.... S m . and assume zero initial conditions, the 
response y{i) is 

>•(/) = y,(f) + y,,(f) 

n m 

= 2>« + £*< 

where r, denotes the residue of pole v, in 
H(.\)X( *)« '■' and 

R, denotes the residue of pole S, in 

If the poles of H(\) and AYvl are all simple 
and distinct, and if the poles of //M do not 
cancel Ihe zeros of .V(v) or vice versa, then all 
residues will be of the form 

r, = r,e'.< and R, = R,e\' 
See the Appendix or reference 2 lor methods 



Of computing residues. 



If the initial conditions are not zero the com- 
plete response can he obtained to within in arbi- 
trary constants. Again assuming simple, distinct 
poles and no pole-zero cancellation, yd) will he 

y(') = i>, +v* t «v. 

The constants A , are determined from Ihe initial 
conditions, which are usually given as values of 
.v(r) and its derivatives at / = 0. 

For many forcing functions and for most 
systems, the forced response y M) is the steady- 
state response. Examples of aperiodic inputs for 
which this is true are step functions, ramp func- 
tions, and suddenly applied periodic functions. 

However, if the time constants of the natural 
frequencies — poles of H(s) — arc greater than 
the time constants of the input frequencies — 
poles of X(s) — then the forced response will 
decay to zero faster than the transient response. 
In some systems, such as those made up of only 
L and C. and no R. the transient response will 
not die out at all. The transient response, there- 
fore, may actually be the total steady-slate 
response. ■ 
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Table II Input-Output Relationships for a 
Linear System. System Function - H(s) 



Type ol Input 


Translorm ot Input 


Transform ol Output 


Output 


Damped Sinusoidal 

X(t) = A C. 1 COS (m,J + </>o) 

In phasor notation: 


X(s) = j*A [e>**($ — s„) 
+ e".B(s — 

4i = — iT.i + /"h> 


Y(s) = H(s)X(.\) 


- Re[HU.,)A efaeV) 


Xt = A 

It is understood that this means 
x(i) ^ Re[A e**:'] 


Y T = H(s„)X, 


yll) = Re[HUJX,e--') 


Periodic 

x(r) = Y,A. cos (n.„„r -f £„) 
In phasor notation: 


Transform of the /ith harmonic 

*„(/,.,) = tA. [<•"».«(«. - t ) 

+- e '«.8( M -f «*)] 


By superposition 

r(M - £ «(/...) *..<>) 


yd) = £««-|«(//» 1. •!.(■"• ..•"-.') 




y. = mjn»,„)x. 


v(r> - Y.RflH(in„,)X„e"^\ 


Aperiodic 

.VI/) 


Y(i...\ f x(t)i- 1 'ill 

t\ YJHII I A\l It III 

■ ' - 

if integral exists, 
or 

X(S)=r J ' X (t)e «dt 


K(/u.) = H(jm)X(J») 
or 

KO) = W(j)AT(s) 


l f " 

yto=^~ y ina,.)X(j,.)w 'dm 

or 


Sampled 

x*(t) = £ x(n)S It - nT) 


XU.) = f; x(n)z • 


where 

Wj) = £ / ""»" " 
where 

li(n) = unit pulse response 


v('D - ~ '/ n>:)X(z)z" <dz 
2-/ J 

C 


Random 

R,(t) — Input Autocorrelation 
Function 


Input Power Spectral Density 
= f R,(r)e '-"dr 


SJfa) — | #(/«.) |* £-(/'..>) 


Output Autocorrelation Function 



Appendix. Complex-Variable Formulas 

1. (?) f(z)dz = 0 if j(z) is analytic on and inside the 
g closed contour C. 

2. n(z)dz = 2-rj ^2 [residues of poles of (Hz) inside C|. 



3. If f(z) is analytic on C, and C, and in the region between 
them, then /(;) can be represented in that region by a 
Laurent Series. 



If y<z) has an nth-order pole at r„ then 
I T d" 

residue of z„ = 



/(;) = J^AJz-Zr,)" 

1 - J 




(»-»! 



efe" 



C is any curve between C, and and enclosing z„. 
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